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NAME:  CHIPING  LI 


BRIEF  DESCRIPTION  OF  PORTFOLIO: 


Meet  Basic  Combustion  Challenges  with  New  Approaches 
Explore  New  Energy  Conversion  Opportunities  for 

Next  Generation  Air  Force  Propulsion  Systems  of  Game-Changing  Efficiency  and  Operability 


Key  Portfolio  Attributes; 

•Understand  Fundamentals  in  Realms  of  Air  Force  Interests  (understand  the  nature  as  it  is) 
•Quantify  Rate-Controlling  Processes  and  Scales  in  Multi-Physics,  Multi-Scale  Phenomena 
(find  ways  to  control  complex  phenomena) 

LIST  SUB-AREAS  IN  PORTFOLIO: 

1.  Combustion  Chemistry  (underlying  chemistry,  new  approaches  --  working  with  Drs.  Berman/RTE,  ARO  &  DOE/BES) 

2.  Turbulent  Flame  Properties/Models  (nonlinear  flow-chemistry  interaction,  new  thinking/tools) 

3.  Combustion  Numerics  (new  tools  --  collaborating  with  Dr.  Fariba/RTA) 

4.  Combustion  Diagnostics  (new  tools  --  collaborating  with  Dr.  Parra/RTB) 

5.  Game-Changing  Energy  Conversion  Concepts  (new  opportunities,  Drs.  Berman/RTE,  Luginsiand/RTB,  &  onr) 
Sub-areas  are  multi-disciplinary:  collaboration  with  other  POs  &  Agencies  are  essential. 
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Combustion  -  the  Central  Process  in  Converting 
...  Chemical  to  Mechanical  Energy  in  AF  Propulsion  Systems 

•  Most  Important  determining  factors  of  operability  and  performance; 

•  Least  Understood  areas  in  basic  combustion  research,  with  large  uncertainties; 

•  Confluence  of  “grand-old”  fundamental  science  challenges,  immediate  needs  and  long-term 
interests. 


Cavity  Based  Scramjet  Combustor 


combustion-Mixing 
interaction 


shock 


Cooling-related  fuel 
Property  issues 


Flameholding  Flame  propagation 
(flammability  limit)  (tu  rbu  lent-f lame  speed) 
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Portfolio  Directions/Trends 


1.  Combustion  Chemistry 

•  Reaction-pathway  centric  approaches  (based  on  ab  initial  methods  and  exp.  data) 

•  Traditional  detailed  reaction-rate-constant  centric  approaches  ^ 


2.  Turbulent  Flame  Properties/Models 

•  Turbulent  flame  experiments  in  realm  of  AF  interests  (high-Re,  compressible) 

•  Laminar  and  weakly  turbulent  flame  experiments 

•  Scale  interaction  models  based  on  DNS,  experimental  data  and  new  math  approach  t 

•  Models  based  on  assumptions  not  directly  verifiable  by  experiments  ^ 

3.  Combustion  related  Numerical  Techniques 

•  Coupled  simulation-experiment  approaches^ 

•  Computational  methods  for  studying  stochastic  pathways^ 

•  Computational  methods  for  reduction  of  “large”  detailed  combustion  chemistry  models^ 


4.  Combustion  Diagnostics— > 


5. 


Game-Changing  Energy  Conversion  Concepts 

•  New  combustion  concepts— ► 

•  Direct/partially  direct  conversion  from  chemical  energy  to  mechanical  energy^ 

•  New  energy-storage/  fuel  concepts  for  propulsion  application  S' 
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Coordination  with  Other  Agencies 


1.  Strong  collaboration  is  continuously  being  forged  in  following  areas: 

•  Diagnostics  (Mainly  DoE,  NASA) 

•  Numerical  (DoE,  NASA,  ARO) 

•  Combustion  Chemistry  (DoE,  ARO,  NSF) 

•  Innovative  Combustion  Concept  (ONR,  ARO) 


2.  Dividing  problems  and  condition  areas  according  to  each  interests: 

•  AFOSR  combustion  portfolio: 

•  Turbulence  combustion  area:  Air-Force  relevant  realms,  i.e.  compressible,  high-Re  conditions  for 
propulsion  applications 

•  Combustion  Chemistry:  Reaction-pathway  centric  approaches 

•  DOE  -  a  well  funded  combustion  program  focusing  on  basic  energy  research: 

•  turbulence  combustion  area:  ground-base  energy  systems  and  auto-engine  types  of  applications  at 
relatively  low-speed  and  low-Re  conditions  (TNF  etc.) 

•  combustion  chemistry:  large,  detailed  reaction-rate-constant  centric  approach 

•  NASA  -  a  modest  combustion  program  focusing: 

•  “Very-high”  speed  (space  access)  region 

•  Overlapping  interests  and  close  coordination  with  AF  programs  (scramjet,  rockets  etc.). 

•  NSF  -  a  modest  combustion  program: 

•  Covers  broad  ranges  of  combustion  problems 

3.  Multi-Agency  Coordinate  Committee  of  Combustion  Research  (MACCCR) 

•  Functioning  well  and  its  positive  roles  will  continue 


Multi-Agency  Collaboration  Benefits  Every  One 
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Combustion  Chemistry 

a  New  Direction 
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Combustion  Chemistry: 
History  and  Recent  Progress 


Spalding  used  3-step  “fictitious” 
reactions  to  describe  the  global 
reaction  kinetics  of  a  flame 
(Philos.  Trans.  Royal  Soc.  London 
A  1956) 

A ►  2 B  —  06  kcal/mole  =  radical-producing, 

B- \-A - ►  B  -f  2C-f  33  kcal/mole  =  chain, 

B+B+X — >2C+X+99  kcal/mole  =  chain-breaking. 


Late  60s-early  70s 

remove  empiricism  by  using  ~10- 
20  step  elementary  chemistry 


burned 

Intermediates 

unburned 


10 


add  physics 
or  just  follow 
the  formula? 


10 


Ranzi  mechanism 
comlete,  ver  1 201 


2-methyl  alkanes  (LLNL)  • 

■ 

methyl  palmitate  (CNRS) 


Biodiesel  (LLNL) 

★ 

Cl 6  (LLNL)' 


Cl  2  (LLNL) 
CIO  (LLNL)-'  I 

ir F 


MD  (LLNL) 

l 

Cl  4  (LLNL) 

PRF  (LLNL) 


^  — —  n-heptane  (LLNL) 
J'etSURF  2.0 


iso-octane  (LLNL) . 

iso-octane  (ENSIC-CNRS)  ^ 

X  T  n-l 

Gasoline  (Raj  et  al) 

CH4  (Konnov)  *’• 

▲  .  B  skeletal  iso-octane  (Lu  &  Law) 

USCC1-C4  B  -M  skeletal  n-heptane  (Lu  &  Law) 
USC  C2H4  •  .  -  '•n-butane  (LLNL) 

■  .t, 3-Butadiene 

j  j  DME  (Curran) 

neo-pentane  (LLNL) 


Cl  -C3  (Qin  et  al) 
GRI3.0  • 


^C2H4  (San  Diego) 
’  CH4  (Leeds) 


I  =  5K 


before  2000 
2000-2004 
2005-2009 
since  201 0 


10 


10 


1950s’ 

Three-step,  global  chemistry 
with  detailed  transport 

■  burned 

unburned 


Graham  Dixon-Lewis 


Dixon-Lewis’ 
calculation  of  the 
structure  of  a 
hydrogen  flame, 
including 
intermediate 
radical 

concentrations, 
using  detailed 
chemistry 
(Proc.  R.  Soc. 
Lond.  >4  1970) 


10  10 
Number  of  species,  K 

Today:  following  the  formula,  an  extrapolation  of  Dixon-Lewis 
work,  but  this  is  sure  not  what  Dixon-Lewis  had  in  mind. 

•  0(1 04)  species  and  0(1 05)  reactions  with  rates  and 
pathways  unverifiable. 

•  Back  to  the  empirical  past:  >  90%  species  considered 
are  not  detectable  by  any  experimental  means. 


Ranked  maximum  concentrations  of  species  computed  as 
an  initial  value  problem  using  a  typical  reaction  mechanism 
(T0  =  1400  K  and  P=  1  atm) 


Wang  USC 
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Number  of  reactions,  I 


Uncertainty  vs.  Model  Size 


1 0  r 


10°  r 


Ranzi  mechanism 
comlete,  ver  1201 


2-methyl  alkanes  (LLNL)  * 

■  • 

methyl  palmitate  (CNRS) 


Biodiesel  (LLNL) 

★ 

Cl  6  (LLNL)', 


■  MD  (LLNL) 

Cl  2  (LLNL)- '  ■ 

C10(LLNLK  ■  Cl  4  (LLNL) 

iso-octane  (LLNL)  .  JF  PRF  (LLNL) 

iso-octane  (ENSIC-CNRS)  ▲  A 

★  1 

Gasoline  (Raj  et  al) 


^  ——  n-heptane  (LLNL) 
JetSURF  2.0 


CH4  (Konnov) 


USC  C1-C4 
USC  C2H4  • 


skeletal  iso-octane  (Lu  &  Law) 
skeletal  n-heptane  (Lu  &  Law) 


•  n-butane  (LLNL) 
t  ,3-Butadiene 


;  Cl  -C3  (Qin  et  al)  _  •  DME  (Curran) 

GRI3.0  #  •  neo-pentane  (LLNL) 

• 

before  2000 

GRI1  2  ^^C2H4  (San  Diego) 

A 

2000-2004 

, 

. '  *  CH4  (Leeds) 

■ 

2005-2009 

10 

1  =  5K 

★ 

...  i 

since  201 0 

io1  io2 

io3 

i  o 


The  larger,  the  better  -  Maybe  Not 

Does  the  model  uncertainty  reduce  as 
the  model  size  grows  -  Not  Necessarily 


•  Uncertainties  of  Model  inputs 

•  Uncertainties  with  model  complexity 

•  Relationship  between  model  size  and 
its  uncertainty 


Find  out  optimal  model  size  for  minimal 
uncertainty  —  new  direction  in  UQ  for 
the  combustion  chemistry  model 


Number  of  species,  K 
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V  ^Two  Distinctive,  Complementary  Approache 

for  Real  HC  Fuel 


Traditional/State  of  Art 
“Rate  Constants  Centric” 

(1)  INCLUDE:  combinatory 
approaches  (up  to  ~  104-105 
reaction  steps  for  common  HC 
fuels) 

(2)  ESTIMATE:  rate  constants 
for  most  and  calculate  and 
measure  for  some 

(3)  SELECT:  through  sensitivity 
analysis,  targeting  ~  order  of  102 
reaction  steps 


Notel :  large  uncertainty  in  steps  1  &  2 

Note2:  There  is  a  confusion  between 
approximate  natures  of  those  large, 
complex  reaction  sets  and  exact 
nature  of  the  human  Genome  set. 


In  the  combustion  process, 
reactions  follow  uncertain 
pathways,  step  by  step. 

k(T)  =  ATe-EalRT 

At  each  step  in  traditional 
Arrhenius  model,  the 
reaction  rate  is  controlled  by 

rate  constants. 


k= A  fexp(-E/RTf 
A  n  E 


references/ 

comments 


Reactions  of  propene 

1  aC3H5  +  H(+M)  =  C3H6(+M) 

2  CH3  +  C2H3  (+M)  =  C3H6  (+M) 


c3h6+h=c2h4+ch3 

C3H6  +  H  =  aC3H5  +  H2 

c3h6+h=ch3cch2+h2 
c3h6 + 0  =  CH2CO  +  ch3  +  H 
c3h6+o  =  c2h5  +  hco 

C3H6  +  0  =  aC3H5  +  0H 

c3h6+o  =  ch3cch2  +  oh 

C3H6  +  0H  =  aC3H5  +  H20 

c3h6+oh  =  ch3cch2+h2o 

C3H6  +  H02  =  aC3H5  +  H202 
C3H6  +  CH3  =  aC3H5  +  CH4 
C3H6  +  CH3  =  CH3CCH2  +  CH4 


2.00xl014 
1.33xl060 
a=0.020  =1097 

2.50xl013 
4.27xl058 
a=  0.175  f**=1341 
1.60xl022 
1.70xl005 
4.00xl005 
1.20xl008 
3.50xl0°7 
1.80xl0u 
6.00xl010 
3.10xl006 
l.lOxlO06 
9.60xl0°3 
2.20x10°° 
8.40X10-01 


-12.0  5968 

/=1097  r*  =6860 

-11.94  9770 

/=60000  /*=10140 
-2.39  11180 

2.5  2490 

9790 
327 


2.5 

1.65 


1.65  -972 
0.7  5880 


0.7 

2.0 

2.0 

2.6 

3.5 

3.5 


7630 

-298 

1450 

13910 

5675 

11660 


kx,d 

h 

e 

KJ 

k0 

e 

1  atm,  g 
[33] 
[33] 
[33] 
[33] 
[33] 
[33] 
[33] 
[33] 
[33] 
[33] 
[33] 


New/Start  Exploring 
“Follow  the  Pathway” 

(1)  SELECT:  identify  important 
pathways  following  PES 

(2)  INCLUDE:  only  include  most 
relevant  ones  -  targeting  no  more 
than  -order  of  102  reaction  steps 

(3)  OBTAIN:  rate  constants  from 
experimental  measure  and  ab  initio 
calculations 


Notel :  understanding  of  initial  fuel  break¬ 
up  is  most  important 

Note2:  made  possible  for  recent  develop 
in  diagnostics  and  ab.  initio  chemistry 
calculation  method 


????  -  Have  Explored 
for  more  than  forty  years 


Combustion  Chemistry  Models  of  limited  reaction  steps  with 
acceptable  uncertainties,  usable  for  reactive  CFD  too] 
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automated 
O  PES/PEL  O 
generation 


PEL-RCCE  algorithm  to  enable  application  of  PES  in 
turbulent  combust^n 

Potential  energy 
treated  as  adaptive 
with  increased 
resolution 

Rate-controlled  ( t  =  ti 
Rdte-controlled  ( t=  t. 
Rate-controlled  ( t  =  t: 

Equilibrium  soln  (t  =  0) 


Chemical  Dynamics 


Ofs 

1 50  fs 

lOOfs 

„  • 

«  * 

2007s 

t*  • 

1 25  fs 

•  ® 

250  fs 

'/  • 

> 

Experimental  Validation 


Time  [jus] 


rssr 

I  I.6C-0S 
|  I.4C-05 
'  l.3e-05 
l.lc-05 
<>.0e-06 
7.2e-06 
5.4e-06 
I  3.6e-06 
|  l.8e-06 

i  0.0c‘00 

Y(CI0H8) 
_  2.8c-04 
I  2.5e-04 
|  2.2e-04 
2.0C-04 
i,  l.7e-04 
l.4«-04 
l.lc-04 
X. 40-05 
S.6e-05 
U  2. Re-05 

I  0.0c '00 


Combustion 
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Diagnostics  in  Combustion  Chemistry 


t=200  fs 


F emtochemistry 

A  +  BC  [ABC]*  ->  AB  +  C 

Zewail,  JPC  A  1 04,  5660  (2000) 


Heart  beat 


Picture  refresh 
Period  of  in  LED  TV 
AM  radio 


Molecular 

rotation 


Computer 
processor  speed 


Electronic 

motion 


Molecular 

vibration 


•wmaamaa# 


10-15 
femtosec 


io-21 

zeptosec 


10-i8 

attosec 


10'9 

nanosec 

io-12 

picosec 

Attosecond 
physics/chemistry 

A  +  B  ^  [AB]*  +  e AB 

Murnane,  Science  322,  1207  (2008) 


I 

IO'6 

microsec 

B 


IO’3 

milisec 


second 


\WoakHHG  .  Strong  HHG 

k  V,  Electron  ♦ 

recombination  > 

'’A— 


NzOi 


N-N  Bond  Length 


m 


Gord  et.  al.,  AFRL 
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With  recent  developments  in  combustion  diagnostics 
(especial  ultra-fast  lased  based  diagnostics)  and  ab.  initio 
chemistry  methods,  we  have  unprecedented  opportunities 
in  combustion  chemistry,  —  AFOSR  is  leading  the  charge 

leading  to  usable  models  with  acceptable  uncertainty  to 
revolutionized  Air  Force  propulsion  system  development. 
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Turbulent  Flame  Property 

in  Air-Force-Relevant  Realms 
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Turbulence  Combustion:  Fundamental  Structures, 
Critical  Scales  and  Relevant  Conditions 


Reaction 

£ones 


Thin  Reaction  Zt 


Broken 
Reaction. 
£on@s  . 


CCTTtfi 


f  Lame!e( 


Thin  Reaction  £one 


La|inJ 

Hamos 


Wrinkled  flamelets 


Corrugated  f  lame  Lets 


Laminar 

flames 


PGC 


uto  Engines 


Hypersonics 

Gas  Turbines 


Rockets 


laminar  flame  wrinkled  flame  ball  Hi-M,  Hi-Re  (flame  brushes) 


(1) Little  Understanding  and  Data  Available  at  AF  Relevant  Compressible,  High-Re  Conditions; 

(2) Needs  for  Better  Definition  of  Re-Conditions  in  Regions  of  Interests 
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High-Re,  Compressible  Turbulence  Flame 
^♦•Experiments  at  AF  Relevant  Condition  Ranges 


•Focus  on  key  combustion  properties  and  characteristics  such  as: 

•Flame  propagation, 

•Flammability  limit 
•Combustion  instability 

•Multi-phase  conditions  applicable  to  Air  Force  propulsion  systems 
•Made  possible  by  diagnostics  developed  by  this  portfolio  up  to  date 


Fuel-air  Final  Mixing 

Mixing  to  uniform, 

Region  homogeneous 

i 

<s*ag7  t 

Grid 

Generated 
Turbulence  to 
Control  (u',  X) 


Pilot  name  / 
plasma 


Key  Requirements  (Experimental  Data  Objectives): 

1.  Understanding  the  above  key  combustion  phenomena  and  characteristics; 

2.  Quantifying  rate-controlling  processes  and  scales  that  govern  those 
phenomena  and  characteristics; 

3.  Developing  and  validating  as  directly  as  possible  basic  model  assumptions 

4.  Controlling  and  quantifying  turbulence  properties  are  essential. 

Proposals  are  being  considered  and  funded  for: 

•Defining  relevant  conditions  and  Studying  Critical  Scales  (1  funded  in  FY12) 
•Relevant  Experiments  in  different  configurations  (4  funded  in  FY12) 


Understanding  Nature  from  Observation  and  Data 


Turbulent  fuel  stream 


r  I V  .  /  , 

-  - 

I  Laser/  *  dy  Flame  * 

PUF  /— ( - - - -  ' 


Turbulentflow  reactor 


Reactive  flow 
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More  Turbulent  Flame  Experiments 


Speed  and  Free  Propagating  Turbulent  Flame  and  Swirl  Flame:  (PI:  Lieuwen,  Georgia  Tech) 


swirl  burner  (jet-engine  similar)  data 


sister  program  Bunsen  burner  data 


o 


CO 


Q 

o 

sT 


co 


- 1 - 

X 

X 

w 

— 1 - 1 - 1 - 1 - 

X 

x  x  X 

+ 

X 

X  ,  X 
x+++ 

X 

X 

X 

_ 1 _ 

+  20  m/s,  CH4,  <t>  =  0.9 

40  m/s,  CH4,  <|>  =  0.9 
x  CH4  (Littlejohn  et  al.  2010) 

_ i _ i _ i _ i _ 

0 


10  12 


u'  /5,n 

ax  L,0 


Major  data  set  obtained  over  range  of 
velocities  (4-70  m/s),  pressures  (1-20  atm), 
turbulence  intensities,  fuel  compositions. 
ST/SLcan  be  »100. 


t  =  0  pis  t  =  1 1 1  m-S  t  =  222  pis  t  =  333  pis  t  =  444  pis 

No  upper  limits  observed  in  turbulence  flame  speed!!! 

Much  more  efficient  and  compacted  combustor  can  be  designed. 
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More  Turbulent  Flame  Experiments 


Flame  Speed  and  Self-similar  Propagation  of 
Turbulent  Premixed  Flames:  (PI:  Law,  Princeton) 


CV:  Check  Valve,  PG:  Pressure  Gauge,  PT:  Pressure  Transducer,  M:  Fan  Motor, 
L:  Cylindrical  Lens,  E:  Electrodes,  W:  Quartz  Window 


tSHr*  ■  W:-  \ 

W  4  -  ‘O 


4  W0-  m 

®  !|£f 


Pressure  =  Pressure  = 

Sstfltleren  images  of  turbulent  premf^JnCH4_air 
flames  ( <f>=0.9 ,  Le=1)  at  same  urms 


Turbulence  flame  speed  can  be  scaled,  at 
least  partially  understood  and  modeled. 
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Turbulent  Flame  Interactions 


Domain  cross-sections 

(a)  Lam.  (b)  FI  (c)  F3  (d)  F5 


Turbulent  Flame  Structure:  (Pis:  Oran,  Poludnenko  &  Hamlington  NRL 


thicken  pre-heating 
(pyrolysis)  zone 


■  small  turbulent  intensity,  1 


large  intensities,  preheat  (pyrolysis)  zone  broadened  and 
reaction  zone  virtually  unaffected 

■  robustness  of  reaction  zone  -  turbulent  diffusion 
suppressed  in  reaction  zone  by  heat  release  (suppressed 
small  scales) 

^  tangential  strain  rate  thins  flame 


We  may  be  able  to  understand  turbulent  flame  structure  after  all... 
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Turbulence-Chemistry  Interaction 


burned 


Flame  zone-  oxidation  of 
cracked  intermediates 

I  Flame  physics  dictates 
decoupled  crackingand 
oxidation 

Preheat  zone  dominated 
by  fuel  cracking 

unburned 


Stage  1 :  quantify  the  pyrolysis  process 
in  the  thicken  pre-heating  zone,  which 
leads  to  ~6  c1-c4  molecular  fragments 

Stage  2:  combined  with  the  c1-c4 
combustion  chemistry  that  has  been 
well  characterized. 


A  partial  inventory  of  available  reaction  model 


10 


10° 


1  0 


Ranzi  mechanism 
comlete,  ver  1 201 


2-methyl  alkanes  (LLNL)  • 

■ 

methyl  palmitate  (CNRS) 


Biodiesel  (LL-NL) 

★ 

Cl  6  (LLNL)', 


1  MD  (LLNL) 

Cl  2  (LLNL)- '  ■ 

CIO  (LLNL)- •  ■  Cl 4  (LLNL) 

iso-octane  (LLNL)  ,  'JT  PRF  (LLNL) 
iso-octane  (ENSIC-CNRS^  AA  n.heptane  (LLNL) 

Gasoline  (Raj  et  al)  JetSURF  2.0 

CH4  (Konnov)  *•' 

A  .  g  skeletal  iso-octane  (Lu  &  Law) 

USCC1-C4  a  . -St  skeletal  n-heptane  (Lu  &  Law) 

USCC2H4  •  .  •  '•n-butane  (LLNL) 

■  .1,3-Butadiene 
A  DME  (Curran) 

neo-pentane  (LLNL) 

~  before  2000 
2000-2004 
2005-2009 
since  2010 


Cl  -C3  (Qin  et  al) 
GRI3.0 


10' 


10° 


10 


Number  of  species,  K 


Turbulence-Chemistry  Interaction:  (PI:  Wang,  USC) 


Turbulent  Pre-Heating  (Pyrolysis)  Zone  Makes  the  Chemistry  Model  Simpler.... 
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With  recent  developments  in  combustion  diagnostics  and 
numerical  simulation  for  the  reactive  flow,  we  begun  to 
observe  and  understand  fundamental  attributes  of  the 
turbulent  flame  in  Air-Force-relevant  realms, 

Leading  to: 

•  Quantify  of  interactions  among  different  scales 

•  Establish  of  usable  turbulent  combustion  models 

•  With  acceptable  uncertainty  to  revolutionized  Air  Force 
propulsion  system  development. 
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Examples  of  Continuous  Transition 


State-of-Art  Optical  Fibers,  Probes,  Single-Beam  Techniques 


Diagnostics 


Basic  Combustion 


J85  Turbojet  Engine 
Equipped 


New  Ignition  Technique 


Shadowgraph  Chemiluminescence 
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Closing  Statements 


After  a  year,  the  portfolio  is  taking  shape. 

Supported  projects  have  started  showing 
very  encouraging  results. 

More  to  come,  stay  tuned. 
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